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ABSTRACT 

The  literature  regarding  transverse  vibrations  of  fluid-conveying  pipes  in  air  is  plenty,  however,  most  of  them 
is  for  the  stability  of  single-span  clamped-free  (C-F)  pipes,  and  the  reports  relating  to  the  critical  flow  velocities  of  the 
multi-span  fluid-conveying  pipes  (either  in  air  or  in  water)  with  various  end  supports  are  not  found  yet,  thus,  this  paper 
tries  to  study  the  title  problem  by  using  the  finite  element  method  (FEM).  First  of  all,  the  property  matrices  of  an 
axial-loaded  fluid-conveying  pipe  element  in  water  are  derived  through  application  of  the  Lagrange’s  equation.  Next, 
the  assembly  technique  is  used  to  construct  the  overall  mass,  damping  and  stiffness  matrices  for  the  entire 
fluid-conveying  pipe,  and  after  imposing  the  total  constrained  degrees  of  freedom  (DOFs),  the  equations  of  motion  for 
the  entire  multi-span  fluid-conveying  pipe  are  obtained.  Solution  of  the  latter  equations  of  motion  will  produce  the 
complex  eigenvalues  and  complex  eigenvectors,  with  the  imaginary  part  of  each  eigenvalue  and  the  real  part  of  each 
eigenvector  denoting  the  damped  natural  frequency  and  damped  mode  shape,  respectively.  To  confirm  the  reliability  of 
the  presented  approach,  the  stability  curve,  the  Argand  diagrams  for  several  complex  eigenvalues  and  the  associated 
dimensionless  critical  flow  velocities  of  a single-span  pipe  with  three  (end)  boundary  conditions  (BCs)  are  determined 
and  compared  with  the  existing  literature,  and  good  agreement  is  achieved.  For  a “q-span”  pipe  ( q = 1,2,3,...  ) with 

equal  span  lengths  and  S end  supports  (or  BCs),  it  has  been  found  that,  its  first  dimensionless  critical  flow  velocity 
tq^U  \ cr  leading  to  buckling  can  easily  be  determined  from  the  Argand  diagram  associated  its  I'1  complex 
eigenvalue  X{,  and  this  is  not  true  for  a classical  axial-loaded  “q-span”  beam  on  the  basis  of  static-stability  theory. 
Furthermore,  for  the  last  pipe  with  pinned-pinned  (P-P)  end  supports,  it  has  been  found  that  the  first  dimensionless 
critical  flow  velocity  is  given  by  ppUlcr  = q7t . In  addition  to  the  stability  analyses,  the  determination  of  damped 

natural  frequencies  and  damped  mode  shapes  of  a multi-span  fluid-conveying  pipe  (with  Coriolis  effect  considered)  in 
three  BCs  is  also  illustrated.  Since  the  information  regarding  the  dynamic  behaviors  of  the  multi-span  fluid-conveying 
pipes  (either  in  air  or  in  water)  is  rare,  the  presented  approach  and  the  obtained  results  will  be  useful  in  the  practical 
applications. 
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1.  INTRODUCTION 

The  dynamic  behaviors  of  the  fluid-conveying  pipes  are  the  important  information  for  the  engineers,  thus,  the 
reports  concerned  are  numerous.  In  the  work  of  Ibrahim  [1],  more  than  four  hundred  reports  regarding  various  fluid- 
conveying  pipes  are  reviewed  and  it  is  found  that  most  of  them  are  for  the  single-span  pipes  with  cantilevered  conditions  as 
one  may  see  from  the  list  of  references.  In  early  years,  Chen  [2],  Paidoussis  and  Issid  [3],  Ting  and  Hosseinipour  [4], 
Paidoussis  et  al.  [5],  Chen  and  Fan  [6],  Pramila  et  al.  [7],  Dupuis  and  Rousselet  [8],  To  and  Healy  [9],  Pramila  [10],  Chang 
and  Cheng  [11],  and  Bratt  [12]  will  be  some  of  the  pioneers  for  the  dynamic  analyses  of  single-span  pipes  and  only  a few 
reports  regarding  the  multi-span  pipes,  such  as  the  works  of  Stein  and  Tobriner  [13],  Singh  and  Mallik  [14,15],  Koo  and 
Park  [16],  Wu  and  Shih  [17]  and  Wu  et  al.  [18].  From  the  existing  literature  one  sees  that  the  dynamic  behaviors  of  most 
single-span  pipes  are  studied  by  using  the  assumed  mode  (or  modal  analysis)  method  [2-4,18]  or  the  FEM  [6-10,18],  but 
those  of  the  periodically  supported  multi-span  pipes  are  investigated  with  the  wave  approach  [14-16],  in  which,  the 
dynamic  equation  for  any  span  of  a periodically  supported  multi-span  pipe  is  assumed  to  be  the  same  so  that  the  CPU  time 
of  computations  is  independent  on  the  total  number  of  spans,  however,  this  assumption  may  disagree  with  the  realistic 
conditions.  In  the  work  of  Wu  and  Shih  [17],  the  transfer  matrix  method  (TMM)  [19],  the  mode  superposition  method  [20] 
and  the  Newmark  direct  integration  method  [21]  are  used  to  determine  the  natural  frequencies  and  mode  shapes  of  the 
undamped  multi-span  pipe  as  well  as  the  forced  vibration  responses  of  the  damped  one  with  Coriolis  force  considered. 
Although  the  formulation  of  [ 17]  is  more  practical,  but  the  stability  analysis  of  the  multi-span  fluid-conveying  pipe  is  not 
performed.  Thus,  in  Ref.  [18],  Wu  et  al.  have  studied  the  last  problem,  but  only  the  3-span  pipe  in  air  and  with  pinned- 
pinned  (P-P)  end  conditions  was  studied. 

From  the  foregoing  literature  review,  it  is  seen  that  the  reports  relating  to  the  critical  flow  velocities  of  the  multi- 
span fluid-conveying  pipes  (either  in  air  or  in  water)  with  “various”  end  supports  (or  BCs)  do  not  appear  yet,  thus,  this 
paper  tries  to  present  some  information  in  this  aspect.  First  of  all,  the  property  matrices  of  an  axial-loaded  fluid-conveying 
pipe  element  in  water  are  derived  through  application  of  the  Lagrange’s  equation  [22-26].  Here,  for  simplicity,  the  main 
effect  of  the  surrounding  water  on  the  fluid-conveying  pipe  is  assumed  to  be  distributed  on  the  added  mass  [27,28].  Next, 
the  assembly  technique  is  used  to  construct  the  overall  mass,  damping  and  stiffness  matrices  for  the  entire  fluid-conveying 
pipe,  and  after  eliminating  all  constrained  degrees  of  freedom  (DOFs),  the  equations  of  motion  for  the  entire  pipe  system  to 

take  the  form  [m]  [if  } + [c]  {if  } + [£]{«}  =0  is  obtained.  For  convenience,  the  last  equations  are  transformed  into  another 

form  [M\{Y } +[K]{Y } = 0 [22]  and  then  the  existing  approach  is  used  to  solve  the  eigenproblem  [29].  Because  of  the 
effects  of  Coriolis  force  and  the  inner  (or  outer)  damping,  the  eigenvalues  and  eigenvectors  of  the  last  equations  are  the 
complex  numbers.  Numerical  examples  reveal  that  the  imaginary  part  of  each  complex  eigenvalue  denotes  the  damped 
natural  frequency  and  the  real  part  of  each  corresponding  complex  eigenvector  denotes  the  associated  damped  mode  shape. 
In  the  stability  analysis,  it  is  found  that,  for  a “(/-span”  fluid-conveying  pipe  ( q = 1,2,3,...  ) with  S end  supports 

( S = P-P,  C-C  (clamped-clamped)  or  C-P  (clamped-pinned)),  if  its  “kth”  buckled  mode  shape  iq)sYk  (x)  takes  the  same  form 
as  the  “rth”  buckled  mode  shape  sZr(x)  of  a classical  axial-loaded  “1-span”  beam  with  the  same  S end  supports  obtained 
from  the  theory  of  static  stability  [30],  i.e.,  <q)sYk  (x)=sZr  (x) , then  the  kth  critical  “dimensionless”  flow  velocity  <q^Ul  of 
the  “g-span”  fluid-conveying  pipe  is  equal  to  the  “rth”  equivalent  one  sVr  of  the  classical  axial-loaded  “1-span”  beam , i.e.. 
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('sUl=sV'.  . For  convenience,  the  statically  applied  axial  load  N = mfV  associated  with  the  rth  critical  dimensionless 

flow  velocity  sVr  of  the  classical  axial-loaded  “1-span”  beam  is  called  the  “equivalent”  axial  load  and  derived  in  the 
Appendix  at  the  end  of  this  paper,  where  mf  denotes  the  fluid  mass  (inside  the  pipe)  per  unit  length  and  V = U denotes 

the  “dimensional”  flow  velocity.  In  addition  to  the  stability  analysis,  the  free  vibration  characteristics  of  a 4-span 
fluid-conveying  pipe  with  three  kinds  of  end  supports  (or  BCs)  are  also  investigated. 

To  confirm  the  presented  approach  and  the  developed  computer  programs  for  this  paper,  the  stability  curve,  the 
Argand  diagrams  of  several  complex  eigenvalues  and  the  associated  critical  “dimensionless”  flow  velocities  of  instability 
(or  buckling)  for  an  “1-span”  fluid-conveying  pipe  with  pinned-pinned  (P-P),  clamped-clamped  (C-C)  and  clamped-free 
(C-F)  end  supports  are  determined  and  compared  with  the  available  literature,  and  good  agreement  is  achieved.  It  is  noted 
that,  the  mode  shapes  corresponding  to  the  “damped”  natural  frequencies  of  either  single-span  or  multi-span 
fluid-conveying  pipe  (with  Coriolis  effect  considered)  are  not  presented  in  most  of  the  existing  literature,  however,  they  do 
appear  in  this  paper. 

2.  FORMULATION  OF  THE  PROBLEM 


This  section  includes  the  assumptions,  the  derivation  of  property  matrices  of  the  axially  loaded  pipe  element  in 
water,  the  constrained  DOFs  for  various  supporting  conditions  and  the  equations  of  motion  for  the  entire  multi-span  fluid- 
conveying  pipe  along  with  their  solutions. 

2.1  Assumptions  for  the  Formulation 

The  formulation  of  this  paper  is  based  on  the  following  assumptions  [2,3]:  (a)  The  pipe  is  made  of  homogeneous 
Kelvin- Voigt  viscoelastic  material;  (b)  The  fluid  is  incompressible  and  inviscid;  (c)  The  effect  of  the  pipe  motion  on  the 
fluid  is  negligible;  (d)  The  velocity  profile  of  the  fluid  is  uniform  over  the  cross-sectional  area  of  the  pipe  (i.e.,  the 
plug-flow  type).  In  addition  to  the  last  assumptions  for  the  fluid-conveying  pipes  in  air,  for  simplicity,  we  also  assumed  that 
the  main  effect  of  the  surrounding  water  on  a fluid-conveying  pipe  is  distributed  on  the  added  mass  [27,28]. 

2.2  Property  Matrices  of  an  Axial-Loaded  Pipe  Element 

Figure  1(a)  shows  the  mathematical  model  of  a g-span  uniform  fluid-conveying  pipe  submersed  in  water  with  P-P 
supports  in  the  global  coordinate  system  xy , in  which,  the  symbols  1,  2,  -.■Si-1.-.  St,...,Si+ 1,-,  n and  n+ 1 dente 

the  numberings  of  nodes,  and  (1),  (2),  ...,  (i),  (i+1),  ...,  (q-1)  and  ( q ) denote  those  of  spans.  Figure  1(b)  shows  a typical 
two-node  4-DOF  pipe  element  subjected  to  the  axial  load  N in  the  local  coordinate  system  xy , in  which,  uy  = uy(x,t) 

denotes  the  element  displacement  at  position  X and  time  t,  u . ( j = 1 to  4)  denote  the  nodal  displacements,  and  l denotes 

the  element  length.  According  to  the  assumptions  given  in  the  last  subsection,  the  kinetic  energy  T and  potential  energy  V 
of  the  fluid-conveying  pipe  element  shown  in  Figure  1(b)  are  given  by  [7,24,26]: 


T = y J (mp  +ma  + mf)u2ydx  + j mfUuyu'ydx  + ^ J mfU1uydx 
V = \^EIu"y1dx-\  |o  ( Pp  +pfAf)u'ydx 


(1) 

(2) 
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Furthermore,  the  non-conservative  virtual  work  done  by  the  damping  forces  is  given  by  [2,3] 

SWm:  =-("  EUi'du’dx-  C u Su  dx  (3) 

Jo  Jo 
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(b) 


Figure  1:  (a)  The  Mathematical  Model  for  a (/-Span  Uniform  Fluid-Conveying  Pipe  in  Water  with 

P-P  Supports  in  the  Global  Coordinate  System  xy , with  1,2, ...,  .S',  Si Si+1 n and  n+1 
Denting  the  Numberings  of  Nodes,  and  (1),  (2), ...,  (t),  (i+1), ...,  (q-1)  and  (q)  Denoting  Those  of 
Spans ; (b)  the  Two-Node  4-DOF  Pipe  Element  Subjected  to  the  Axial  Load  N in  the  Local 
Coordinate  System  xy,with  u , =uy(x,t)  Denoting  the  Element  Displacement  and 

u j ( j = 1 to  4)  Denoting  the  Nodal  Displacements 


In  Eq.  (1),  mp,  ma  and  mf  are  the  pipe  mass,  added  mass  and  fluid  mass  per  unit  length,  respectively,  U is  the 
constant  axially  flow  velocity,  u'y  and  uy  are  the  differentiations  of  ujx.t)  with  respect  to  coordinate  x and  time  t, 

respectively.  In  Eq.  (2),  E and  I are  Young’s  modulus  and  moment  of  inertia  of  cross-sectional  area  of  the  pipe  wall, 
respectively,  P is  the  externally  applied  axial  force  on  the  pipe  wall  with  the  compressive  one  to  be  positive  (+),  pf  is  the 

pressure  intensity  of  fluid  (inside  the  pipe),  and  A,  is  the  inner  cross-sectional  area  of  pipe  (or  cross-sectional  area  for 
flowing  fluid).  Note  that,  in  Eq.  (1),  the  2nd  and  3ld  terms  are  due  to  the  Coriolis  force  and  centrifugal  force  of  fluid, 
respectively.  In  Eq.  (3),  E is  the  coefficient  of  internal  dissipation  of  the  pipe  material  and  CD  is  the  coefficient  of  external 
viscous  damping  due  to  friction  between  the  pipe  and  the  flowing  fluid. 

According  to  Refs.  [20,23],  one  has  the  following  relationships 

uy  =[ay][u]  (4) 

where 


{u}  = {u1  u2  u3  m4}  = [«,  u2  u3  n4]r,  [ay]  = [aui  au2  au3  au4]  (5a,b) 

In  Eq.  (5b),  [ a ] is  the  shape  function  matrix  with  its  coefficients  given  by 
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aul{%)  = l-3f +2<f , aa2{&  = (t- au3(&  = 3<f-2<f , fl.4(f)  = (-<f  + <f ) t (6a-d) 

with 

£ = x/t  (7) 

Where  x is  the  axial  local  coordinate  and  l is  the  element  length  as  shown  in  Figure  1(b). 

From  Eqs.  (3)  and  (4),  one  obtains  the  non-conservative  generalized  force  vector 

{Qnc } = -f°EI[a;]T[a;]dx{,i}-j°C0[ay]T[ay]dxW}  (8) 


Now,  substituting  Eq.  (4)  into  Eqs.  (1)  and  (2),  respectively,  then  introducing  the  resulting  expressions  for  kinetic 
energy  T and  potential  energy  V as  well  as  the  generalized  force  { Q„,  1 given  by  Eq.  (8)  into  the  next  Lagrange’s  equation 
[22]  with  z = Zj  ={u] 


d_ 

dt 


dT  | dT  dV  ... 

-3—^ +3—  = Qr,nc’  r = l,2,3,...,n 

oz.  j dzr  oz. 


One  obtains  the  equations  of  motion  for  the  pipe  element  shown  in  Figure  1(b)  to  be 
[m\e{u}e+([Cct  ~[Cc]eme  +[*]>}.  =-[Cd]e{d)e 


or 


[m]e{u}e+[c]e{u}e+[k]e{u}e  = 0 
Where  the  subscript  “e”  refers  to  the  pipe  “element”  and 

re  T 

['«]„  =(inp+ma+mf)j[ay]  [ay]dx 


[c\„  = ([Cj;,  ~[Cc]e)  + [Cd]e 

[k]e=[kB]e-[kG]e 

with 


(9) 

(10) 

(11) 

(12) 

03) 

(14) 


[Cc]e=mfU^[a'y]T[ay]dx,  [C d]e  = b^Wyf [ay]dx  + C0^\ay]T [ay]dx  (15a, b) 

[kB]=EI^[a;]T[a;]dx  , [kcl  = N^[aJ[ay]dx , N = Pp+  mfU2  + pfAf  (16a-c) 

In  Eq.  (16c),  N denotes  the  total  effective  axial  load  due  to  the  externally  applied  force  P , the  centrifugal  force 
mfU 2 and  the  fluid  pressure  p , A , . For  the  case  of  Pp  = p , A , = 0 , the  buckling  of  the  fluid-conveying  pipe  is  due  to  the 
centrifugal  force  m ,U2  and  the  corresponding  value  of  U is  called  the  critical  flow  velocity  Ucr. 
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2.3  Constrained  DOFs  for  Various  BCs 

For  the  multi-span  fluid-conveying  pipe  composed  of  n elements  and  h + 1 nodes  as  shown  in  Figure  1(a),  its 
constrained  DOFs  associated  with  the  BCs  are  as  follows: 

• For  the  pinned-pinned  (P-P)  pipe,  one  has  u]  = 0 at  the  first  node  1 and  w2n+1  = 0 at  the  final  node  n + 1 . Thus, 
the  numberings  for  the  constrained  DOFs  are:  1 and2«  + 1 . 

• For  the  clamped-clamped  (C-C)  pipe,  one  has  = (7,  = 0 at  the  first  node  1 and  u2n+l  = u2n+2  = 0 at  the  final 
node  77  + 1 . Thus,  the  numberings  for  the  constrained  DOFs  are:  1,2,  2?;  + 1 and  2«  + 2 . 

• For  the  clamped-pinned  (C-P)  pipe,  one  has  ul=u1=Q  at  the  first  node  1 and  u2n+1  = 0 at  the  final  node  n + 1 . 
Thus,  the  numberings  for  the  constrained  DOFs  are:  1,2  and  2«  + 1 . 

• For  the  clamped-free  (C-F)  pipe,  one  has  ul=u2=  0 at  the  first  node  1 only.  Thus,  the  numberings  for  the 
constrained  DOFs  are:  1 and  2. 

In  addition  to  the  foregoing  constrained  DOFs  associated  with  the  end  supports  (or  BCs)  of  the  entire  pipe,  there 
is  one  DOF  at  each  intermediate  pinned  support  i is  also  constrained.  If  the  numbering  for  the  above  pinned  support  i is 
S{  as  shown  in  Figure  1(a),  then  the  numbering  for  the  corresponding  constrained  DOF  is  “ 25)  —1  ”.  It  is  evident  that,  for  a 
“g-span”  pipe,  there  are  q - 1 intermediate  pinned  supports  and  the  q - I DOFs  corresponding  these  supports  must  also  be 
eliminated  for  obtaining  the  effective  property  matrices  of  the  entire  pipe  system,  in  addition  to  the  DOFs  corresponding  to 
each  specified  BCs. 

2.4  Equations  of  Motion  for  the  Entire  Pipe  System  and  Their  Solutions 

Assembling  all  the  element  property  matrices  obtained  from  previous  Subsection  2.2,  one  can  determine  the 
overall  mass  matrix  [m],  damping  matrix  [c]  and  stiffness  matrix  [k\  for  the  entire  multi-span  fluid-conveying  pipe,  and 
eliminating  the  rows  and  columns  associated  with  the  constrained  DOFs  for  the  specified  BCs  and  the  intermediate  pinned 
supports  shown  in  the  last  Subsection  2.3,  one  may  obtain  the  system  equations  of  motion  to  take  the  form 

[/7F]{m}  + [c]{m}  + [C]{m}  = 0 (17) 

where,  [in] , [c]  and  [k]  are  the  overall  effective  mass,  damping  and  stiffness  matrices  of  the  entire  pipe  system, 
obtained  from  [m],  [c]  and  [k]  by  eliminating  the  rows  and  columns  associated  with  all  the  constrained  DOFs,  respectively. 

For  free  vibrations,  one  has 

{«}  = {«(0}  = {lV"  (18) 

Where  { Y } is  the  amplitude  vector  of  {u(t)} , a is  the  “damped”  natural  frequency  of  the  entire  pipe  system,  t is 

time  and  j = 4—1  . In  order  to  use  the  existing  approach  [29]  to  solve  Eq.  (17),  one  requires  to  transform  the  latter  into  the 
form  below  [22] 

[M]{F}  + L/n{T}  = 0 (19) 
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Where 


[M] 


[0] 

[m] 


-['«] 

[0] 


(20a-d) 


With  respect  to  each  specified  fluid  velocity  U , the  rth  eigenvalue  cor  and  the  associated  rth  eigenvector  [Yr } of 
Eq.  (17)  take  the  pairs  of  conjugate  complex  numbers 


tor=wriR±jarJ,  {Yr}  = 


(21a,b) 


The  imaginary  part  of  cor , (0rI , denotes  the  rth  natural  frequency  and  the  corresponding  real  part  of  j K } , { Yr  R } , 
denotes  the  rth  mode  shape. 

It  is  noted  that,  for  a single-span  cantilevered  (C-F)  pipe  with  P + p fAf  = 0,  the  following  momentum  transport 
at  the  downstream  “free  end”  must  be  considered:  mfU~uy  +mfUuy  lr=L  [6,24-26].  In  other  words,  for  a cantilevered  pipe 
consisting  of  n elements  and  n + 1 nodes  (cf.  Figure  1(a)),  its  total  DOFs  is  2(n  + l)  before  imposing  the  constrained 
DOFs,  thus,  one  must  add  k2n+12n+2  =mfU1  and  c2n+19„+1  =mfU  to  the  corresponding  coefficients  of  the  overall  stiffness 

matrix  | k | and  overall  damping  matrix  | cj , respectively.  It  is  evident  that  the  effects  of  mfU~  and  m,U  on  the  pinned  or 
clamped  BC  at  right  end  of  the  pipe  are  nil,  if  the  (2n+l)th  DOF  is  constrained. 

The  relationship  between  the  rth  dimensionless  frequency  kr  and  the  corresponding  dimensional  one  cor  is  given 
by 

Ar  =corL2yj  mtot  / (£7)  (22a) 


Where 


m,0,  =mp+ma+mf 


(22b) 


And  for  convenience,  the  following  dimensionless  parameters  are  introduced 


m'  = — ^ -,  IT  = ULJ — — , E'  = 
El 


1 r * _ CM 

' 9 ^n~ 


4m,o,EI 


(23a,b,c,d) 


In  Eq.  (22b),  the  symbol  mtot  denotes  the  total  mass  per  unit  length.  It  is  composed  of  three  components  given  by 


mP  = yPpTtidl  -d?) , ma  =\Cap„ril’  mf  =jPf (24a, b,c) 

Where  pp  , pir  and  p f are  mass  densities  of  pipe  (material),  water  (outside  the  pipe)  and  fluid  (inside  the  pipe), 
respectively;  da  and  dt  are  the  outer  and  inner  diameters  of  the  pipe,  respectively;  and  Ca  is  the  added  mass  coefficient. 
For  the  circular  cross-section,  one  may  set  Ca  =1.0  [27,28]. 
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3.  NUMERICAL  RESULTS  AND  DISCUSSIONS 

For  simplicity,  we  set  E = C0  = 0 in  this  paper,  and  before  the  stability  and  free  vibration  analyses  of  the 
multi-span  fluid-conveying  pipes,  some  numerical  results  are  compared  with  the  available  existing  literature  first. 
The  dimensions  and  physical  properties  of  the  pipes  are:  outer  diameter  d0  = 0.254  m,  inner  diameter  di  = 0.2286  m, 

Young’s  modulus  E = 2.068xl0n  N/m2  and  mass  density  of  pipe  pp  =7850  kg/m3  . The  mass  density  of  water  is 

assumed  to  be  pw  =1000  kg/m3  and  the  mass  density  of  fluid,/?,  , corresponding  to  each  specified  mass  ratio  m'f  is 
determined  from  Eqs.  (22b),  (23a)  and  (24a, b,c).  Note  that  the  fluid-conveying  pipes  are  in  water  unless  it  is  indicated  that 
“in  air  (with  pw  = 0 )”. 

3.1  Validation  of  Presented  Method  and  Developed  Computer  Program 

3.1.1  Free  Vibration  of  the  Axial-Loaded  Multi-Span  Pipe  Filled  with  Still  Fluid 

For  a multi-span  pipe  with  equal  span  lengths  i0  and  hinged  end  supports  in  air  (with  p„.  =0),  and  filled  with 
still  fluid,  its  lowest  buckling  load  Px  cr  and  fundamental  frequency  a>0  , are  given  by  [31] 

Plcr=(x/(Q)2EI,  to0l=(7r/(0)2^EI/(ppAp+PfAf ) (25a, b) 

If  pf  =1000  kg/m3  and  f0=6m,  then  Px  cr=  3.983682xl06  N and  <jj01  = 96.77358  rad/sec.  For  the  case  of 
U = 0 and  pf  = 0,  Table  1 shows  the  lowest  four  natural  frequencies  of  the  P-P  pipe  with  1-,  2-,  3-  and  4-  span  and 
subjected  to  the  axial  load  N = 0.6  P,  c = 2.3902092xl0sN,  respectively,  as  shown  in  Figure  2(a)-(d).  In  Table  1,  coN  r and 
w0  r denote  the  rth  natural  frequencies  of  the  pipe  with  and  without  subjecting  to  the  axial  load  N,  furthermore, 
a>N  r /q>o  r denotes  the  frequency  ratio  between  coN  r and  o)0r  obtained  from  the  presented  FEM  and  [31],  respectively. 
Among  the  values  of  a>Nr/co0r  obtained  from  [31],  those  with  four  decimal  places  (e.g.,  coNr/o)0r  = 0.8436)  are  given  by 
Tables  1-5  of  [31],  and  those  with  two  decimal  places  (e.g.,  0)Nr/a)0r  = 0.6300  =0.63)  are  given  by  Figure  5-8  of  [31]. 
From  the  Table  1 one  sees  that  the  frequency  ratios  ( a>N  r /a>0  r ) obtained  from  the  presented  FEM  (with  number  of  pipe 
elements  “for  each  span”  to  be  nx  = 50  ) are  in  good  agreement  with  the  corresponding  ones  obtained  from  [31]. 


Figure  2:  The  Four  Uniform  P-P  Pipes  (in  air)  Filled  with  Still  Fluid  and  Subjected  to  Axial  Load  N with  Equal 
Span  Lengths  = 1 2 = i 3 = i 4 = £0 : (a)  1-span  with  L = 1 0 ; (b)  2-Span  with  L = 2t 0 ; (c)  3-Span  with 
L = 3f  0 ; (d)  4-Span  with  L = 4C0.  Where  L is  the  Total  Length  of  the  Entire  (Multi-Span)  Pipe 
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Table  1:  Influence  of  Axial  Load  N = 0.6  Pl  cr=  2.3902092xl06  N on  the  Lowest  Four  Natural  Frequencies  (rad/sec) 
of  the  Entire  pipe  (in  Air)  Filled  with  Still  Fluid  and  with  1 to  4 Spans  in  the  Hinged  Supporting 
Conditions  as  Shown  in  Figs.  2(a)-(d)  (with  U = 0 , pf  = 0 , p„  = 0 , p f = 1000  kg/m3  and  i0=  6 m) 


No.  of 
spans 

9 

Total 
length 
L (m) 

Methods 

Frequency 
(or  ratio) 

No.  of  pipe 
elements 

ne  = ni  q 

Mode  No.,  r 

1 

2 

3 

4 

1 

Figure2(a) 

6 

FEM 

X, 

50x1 

96.7736 

387.0944 

870.9630 

1548.3817 

X, 

61.2050 

356.8834 

841.4303 

1519.0723 

FEM 

X/X 

0.6325 

0.9220 

0.9661 

0.9811 

Ref.  [31] 

0.6300 

0.9200 

0.9700 

— 

2 

Figure2(b) 

6x2 

FEM 

®o. 

50x2 

96.7736 

151.1788 

387.0944 

489.9163 

X 

61.2050 

127.5324 

356.8834 

464.3150 

FEM 

0>N,rl  <»0,r 

0.6325 

0.8436 

0.9220 

0.9477 

Ref.  [31] 

0.6325 

0.8436 

0.9220 

— 

3 

Figure2(c) 

6x3 

FEM 

®0,r 

50x3 

96.7736 

124.0168 

181.0902 

387.0944 

(0N, 

61.2050 

96.3617 

160.4495 

356.8834 

FEM 

(0N,/(Oo, 

0.6325 

0.7770 

0.8860 

0.9220 

Ref.  [31] 

0.6325 

0.7700 

0.8800 

— 

4 

Figure2(d) 

6x4 

FEM 

50x4 

96.7736 

112.8974 

151.1788 

195.3323 

(On, 

61.2050 

82.7747 

127.5324 

176.0205 

FEM 

(On,  1(0 „ 

0.6325 

0.7332 

0.8436 

0.9011 

Ref.  [31] 

0.6325 

0.7300 

0.8400 

— 

d a>„  r = rth  natural  frequency  of  the  pipe  without  axial  load  (i.e.  N = 0 ), 


b coN  r = rth  natural  frequency  of  the  pipe  with  axial  load  N = 0.6P,  cr . 

3.1.2  Stability  Analysis  of  the  1-Span  C-F  Pipe 

The  dimensions  and  physical  properties  of  this  1-span  cantilevered  pipe  are  the  same  as  those  for  Table  1 and  its 
stability  curve  is  shown  in  Figure  3(a)  with  ordinate  denoting  the  dimensionless  critical  flow  velocity  U cr  and  abscissa 

denoting  the  mass  ratio  mf  . Besides,  Figure  3(b)  shows  the  relationship  between  dimensionless  critical  frequency  (Ar , )cr 

and  mass  ratio  mf . In  the  two  figures,  the  solid  curves  ( ) are  obtained  from  the  presented  FEM,  and  the  dotted  curves 

( ) are  obtained  from  [32]  for  the  C-F  fluid-conveying  pipe  in  air  (with  pw  =0).  Good  agreement  between  the  solid 

curves  and  dotted  curves  confirms  the  availability  of  the  presented  FEM  and  the  developed  computer  program  for  this 
paper. 
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Figure  3:  The  Relationships  between  (a)  The  Dimensionless  Critical  Flow  Velocity  \J*cr  and 

Mass  Ratio  mf , and  (b)  The  Dimensionless  Critical  Frequency  (Ar , ) „ and  Mass  Ratio 
nij , for  the  1-Span  Cantilevered  Fluid-Conveying  Pipe  in  Air 


3.2  Stability  Analyses  of  the  1-  to  4-  Span  Pipes 

For  convenience,  in  the  subsequent  subsections,  the  span  lengths  1 1 ( / = 1 to  q ) for  each  multi-span  pipe  are 
assumed  to  be  equal  to  each  other  with  (j- = <?D  = 6.0  m,  and  the  flow  velocity  and  eigenfrequency  refer  to  the 
“dimensionless”  ones,  respectively,  unless  particularly  mentioned.  The  Argand  diagrams  for  the  1-  to  4-  span  pipes 
(cf.  Figure  2(a)-(d))  with  three  BCs  are  studied  in  this  subsection:  P-P,  C-C  and  C-P.  According  to  Timoshenko  and  Gere 
[30],  a beam  will  be  not  operated  normally  after  the  onset  of  its  first-mode  buckling,  thus,  in  practice,  the  information 

regarding  the  first-mode  buckling  is  most  significant.  However,  for  comparisons,  the  higher  critical  flow  velocity  Urcr 
associated  with  the  rth  complex  eigenvalues  ( r = 1,2 ) of  each  fluid-conveying  pipe  is  also  studied  in  addition  to  the  first 
one  U \cr . It  is  noted  that,  in  the  Argand  diagrams  shown  in  Figure  4 (for  P-P  pipes),  7 (for  C-C  pipes)  and  10  (for  C-P 
pipes),  the  digit  nearby  each  symbol  • (regarding  first  eigenvalue)  or  X (regarding  second  eigenvalue)  denotes  the 
associated  value  of  flow  velocity  U and  the  spacing  between  any  two  adjacent  symbol  points  is  AU"  =0.2  except  for  the 
first  spacing  to  be  AU  =0.05  (from  U =0.05  to  0.1).  Besides,  some  symbols  are  introduced:  (i)  The  rth  critical  flow 
velocity  of  a fluid-conveying  pipe  (with  P = pfAf  =0  and  JJ  0)  is  denoted  by  (|f/rcr  and  the  associated  “buckled” 
mode  shape  by  <'q)sYr  (x)  with  the  left  superscript  “q"  denoting  the  total  number  of  spans  for  the  entire  pipe  and  the  left 
subscript  “S”  denoting  the  supporting  conditions.  For  example,  pplj  (x)  denotes  the  1st  buckled  mode  shape  of  a 3-span  P-P 
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pipe.  If  only  one  critical  flow  velocity  is  associated  with  the  first  (or  second)  complex  eigenvalues,  then  it  is  denoted  by 
[/,  cr  (or  U2  cr ).  However,  if  the  total  number  of  critical  flow  velocities  associated  with  the  first  (or  second)  complex 

eigenvalues  is  greater  than  one,  then  they  are  denoted  by  U]a  cr , Ulh cr , Ulc cr , ...  for  those  associated  with  the  1st  complex 
eigenvalues  and  by  U2a  cr , Ulhcr , U2ccr,  ...  for  those  associated  with  the  2nd  complex  eigenvalues,  (ii)  The  “equivalent” 
rth  critical  flow  velocity  for  a classical  axial -loaded  “1-span”  beam  (with  axial  load  N = mfV 2 ) obtained  from  theory  of 

static  stability  [30]  is  denoted  by  sVrcr  and  the  corresponding  “buckled”  mode  shape  is  denoted  by  sZrcr(x)  as  one  may 
see  from  Eqs.  (A.4),  (A.5),  (A.10),  (A.ll),  (A.16)  and  (A. 17)  in  the  Appendix  at  the  end  of  this  paper.  It  is  noted  that 
1/  and  Y (obtained  from  dynamic  stability)  are  replaced  by  V and  Z (obtained  from  static  stability),  respectively,  to 
avoid  confusion.  Furthermore,  the  left  superscript  “( q )”  is  also  removed  from  V and  Z , because,  based  on  the  theory  of 
static  stability  for  the  buckling  of  a classical  axial-loaded  beam  [30],  only  the  values  of  V and  Z for  “1-span”  beams 
(q  = 1 ) can  be  obtained. 

3.2.1  With  P-P  end  Supports  (or  BCs) 

The  coordinates  ( AIR  , A,,)  and  ( A,  R , A.  r ) for  some  relevant  points  on  the  loci  of  the  1st  and  2nd  complex 

eigenvalues  of  the  1-span  fluid-conveying  P-P  pipe  are  listed  in  Table  2,  and  the  associated  Argand  diagrams  for  the 
1-,  2-,  3-  and  4-  span  ones  are  shown  in  Figure  4 (a-1)  and  (a-2),  (b-1)  and  (b-2),  (c-1)  and  (c-2),  and  (d-1)  and  (d-2), 
respectively.  In  the  latter  figures,  “Figure  4(q-r)”  denotes  the  Argand  diagram  of  the  ‘//-span”  pipe  associated  with  the 
“rth”  complex  eigenvalues  for  various  flow  velocities  lpJU*r , with  span  number  q = a = 1 , q = b = 2 , q = c = 3 , q = d = 4 
and  r = 1,2  . For  example,  “Figure  4(b-l)”  denotes  the  Argand  diagram  for  the  2-span  pipe  associated  with  the  1st  complex 
eigenvalues  for  flow  velocities  'fpUl  =0.05-11.0.  From  Table  2 and  Figure  4(a-l)  one  sees  that  there  are  three  critical 
flow  velocities  to  be  associated  with  the  1st  complex  eigenvalues  A, : ppUla  cr  = 3.1416  ~ 7t , ppUlbcr  = 6.2831=  An  and 
PpU*lCcr  = 9.4 1 1 8 ~ 3^r . If  the  flow  velocity  ppU[  is  greater  than  the  1st  or  the  3rd  critical  ones,  i.e.,  ppUi>ppU\acr  or 
PpUl>ppUlccr , the  1st  frequency  A, ; will  be  equal  to  zero  (i.e.,  A, , = 0)  and  the  pipe  will  be  buckled,  however,  for  the 
flow  velocity  ppU]  to  be  greater  than  the  2nd  critical  one  ppUlbcr  and  less  than  the  3rd  critical  one  ppU*lccr  , i.e., 

PpU^b  Cr<ri'U>\  < ppUiccr*  Table  2 and  Figure  4(a-l)  and  (a-2)  reveal  that  Au  = A^ , and  AIR  =-A1R,  this  result  indicates 
that,  in  the  last  velocity  range,  the  pipe  will  perform  first-and-second-mode-coupled  oscillations  after  it  restores  stability 
from  its  buckled  state  at  the  2nd  critical  flow  velocity  ppUlbcr.  For  convenience,  the  flow  velocities  ppU*  (r  = 1,2) 

associated  with  the  1st  and  2nd  complex  eigenvalues  in  the  last  velocity  range,  ppUlbcr<ppU * < ppUlccr,  are  denoted  by 

ppUlcoupied  and  ppUi, coupled  > respectively. 

In  other  words,  according  to  Table  2 and  Figure  4 (a-1)  and  (a-2),  the  value  of  At , (i.e.,  the  first  natural  frequency 
of  the  pipe)  decreases  with  the  increase  of  the  flow  velocity  p,',U~  and  approaches  zero  when  = n , which  is  the  first 
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critical  flow  velocity  denoted  by  ppUla  rr  for  buckling  of  the  1-span  P-P  pipe.  After  that,  the  values  of  2, , are  equal  to  zero 
until  ppU[  ~ 2 n , which  is  the  second  critical  flow  velocity  ppU lbcr,  and  for  the  value  of  ppU * greater  than  2 n and  less 
than  3 K , i.e.,  7.K<ppU\  < 3tt  , the  locus  of  the  first  eigenvalues  (denoted  by  symbols  •)  and  that  of  the  second 
eigenvalues  (denoted  by  symbols  X)  combine  on  the  vertical  imaginary-axis  (for  2, , and  22 , ) when  ppU[  ~ 6.3  and  then 
the  loci  leave  the  imaginary-axis  at  symmetric  points,  this  indicates  the  onset  of  coupled-mode  flutter  [3]. 


Table  2:  The  Coordinates  ( 2, . R , Ar , ) for  Some  Relevant  Points  on  the  Loci  of  the  1st  and  2nd  Complex  Eigen 
Values  (for  Various  Flow  Velocities  ppUr  = 0.05  — 1 1 .0 ) of  the  1-Span  Fluid-Conveying  P-P  Pipe 


PPU  r 

0.05 

3.1415 

3.1416 

WU' 

PP  ^ 1 a,cr 

... 

6.2831 

mU‘ 

PP  ^ 1 b,cr 

6.2832 

6.4 

... 

9.4117 

9.4118 

PP  ^ lc,cr 

9.6 

... 

11.0 

A 

^l.R 

0 

0 

0 

0.048 

0.184 

0 

6.411 

10.876 

1.141 

-1.554 

-1.885 

K, 

9.868 

0.072 

0 

0 

0 

0.116 

8.218 

0.231 

0 

0 

0 

0 

22 

^ l.R 

0 

0 

0 

0 

0 

0 

0 

-6.411 

-10.876 

-1.105 

-9.655 

28.652 

'k/ 

39.477 

34.800 

34.799 

9.912 

9.905 

8.218 

0.231 

0 

25.028 

29.832 

(a  - 1) 


Real  part  of  eigenvalues,  >.2H 

(a -2) 


Figure  4:  The  Argand  Diagrams  for  the  1st  and  2nd  Complex  Eigen  values  of  (a)  1-,  (b)  2-,  (c)  3-  and  (d) 

4-  Span  Fluid-Conveying  “P-P”  pipes  shown  in  Figure  2 (a)-(d)  with  Mass  Ratio  m f = 0.5  . In  other  words, 

(a-1)  and  (a-2)  are  Those  of  1-Span  Pipe;  (b-1)  and  (b-2)  are  Those  of  2-Span  Pipe;  (c-1)  and  (c-2)  are 
Those  of  3-Span  Pipe;  and  (d-1)  and  (d-2)  are  those  of  4-Span  Pipe 
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(b-1)  (b-2) 


(c-1) 


Real  part  of  eigenvalues,  ).2JK 

(c-2) 


Figure  4:  (Continued) 
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Real  part  of  eigen  values, 

(d-1) 


Real  part  of  eigenvalues,  ^ r 

(d  -2) 


Figure  4:  (Continued) 
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(e) 
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me  “6  4 
«•«* 


— Coupieo  moaesnapeat  £\j.~a 
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Ax  a :oora rates  x/l_ 


Figure  5:  The  first  to  Third  Critical  Buckled  Mode  Shapes  for  the  (a)  1-,  (b)  2-,  (c)  3-  and  (d)  4-  Span 

Fluid-Conveying  “P-P”  Pipes  with  Mass  ratio  m*f  = 0.5 , and  (e)  the  “coupled”  Mode  Shapes  of  the 

1-span  P-P  Pipe  at  Flow  Velocity  ?pUlcovled=$U*2coupled  = 6.4 . 
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It  is  noted  that,  on  the  horizontal  real  ( Al  R ) axis  of  Figure  4(a-l),  the  space  from  AIR  ~ -0.457  to  -1.885  is  very 
small  so  that  many  symbols  (•)  in  that  space  overlap  each  other. 

The  Argand  diagrams  associated  with  the  1st  and  2nd  eigenvalues  of  the  2-,  3-  and  4-  span  P-P  pipes  shown  in 
Figure  4(b-l)  and  (b-2),  (c-1)  and  (c-2),  and  (d-1)  and  (d-2)  are  much  simpler  than  those  of  the  1-span  pipe  shown  in 
Figure  4(a-l)  and  (a-2).  From  Figure  4(b-l)-(d-2)  one  sees  that  the  first-mode  motion  is  not  coupled  with  the  second-mode 

motion,  and  the  values  of  A1I  decrease  with  the  increase  of  ppt/j*  and  equal  to  zero  when  (pp,U]  >qK  with  q = 2,3,4  . 
Thus,  the  first  critical  flow  velocities  for  the  2-,  3-  and  4-  span  fluid-conveying  P-P  pipes  are  2k  , 3 K and  4 n , 
respectively.  This  is  a reasonable  result,  because  the  “first”  critical  buckled  mode  shapes  ^Y^)  (denoted  by  the  curves 
— • — , with  % = x/L)  of  the  2-,  3-  and  4-  span  fluid-conveying  pipes  (based  on  ppUlcr  =qK  and  A{1  = 0)  shown  in 

Figure  5 (b),  (c)  and  (d),  respectively,  take  the  same  forms  as  the  second,  third  and  fourth  buckled  mode  shapes  PPZrcr( £) 
(with  r = 2,3,4)  of  the  classical  axial-loaded  “1-span”  P-P  beam  shown  in  Figure  A. 1(a)  (in  Appendix)  with  its  rth 
equivalent  critical  flow  velocity  of  buckling  given  by  PPVrcr  = rK  (with  r = 2,3,4)  as  one  may  see  from  Eq.  (A.4) 
(in  Appendix). 

It  is  noted  that  the  two  coupled  mode  shapes  for  the  1-span  P-P  pipe  associated  with  the  same  flow  velocities 
ppU \ ,:,mpi,d~pr  ^2  cmw/,', i =6.4  shown  in  Figure  5(e)  are  not  identical,  because  a mode  shape  is  dependent  on  its  eigenvalue 
Ar  and,  at  ppU lcoupled  = ppU2coupled  = 6.4,  the  real  part  of  the  1st  eigenvalue  A,  is  “opposite”  to  that  of  the  2nd  eigenvalue 
A7 , i.e.,  A,  R=  6.411  = - A 2 R , in  spite  of  the  fact  that  the  imaginary  parts  of  the  last  two  eigenvalues  are  equal  to  each  other 
(i.e.,  A, , = A2I  = 8.218  ) as  one  may  see  from  column  10  of  Table  2. 

3.2.2  With  C-C  end  supports  (orBCs) 

For  the  C-C  BCs,  the  configurations  of  the  1-  to  4-  span  fluid-conveying  pipes  with  equal  span  lengths  ic  =6m 
are  shown  in  Figure  6 and  the  Argand  diagrams  associated  with  the  1st  and  2nd  complex  eigenvalues  of  the  “f/-span”  pipe 
are  shown  in  Figure  l{q-V)  and  (q- 2),  respectively,  with  span  number  q = a = 1,  q = b=  2,  q = c = 3 and  q = d = 4. 
The  coordinates  ( A,  R , Al , ) and  (A2R  , A2I ) for  some  relevant  points  on  the  loci  of  the  1st  and  2nd  complex  eigenvalues  of 
the  1-span  and  3-span  C-C  pipe  are  listed  in  Tables  3(a)  and  (b),  respectively.  From  Table  3(a)  and  Figure  7(a-l)  one  sees 
that  all  the  three  critical  flow  velocities,  ccU\acr  =6.2832=  2k  , cc^lbcr  =8.9868  and  ccU\c  cr  = 1 2.5665  ~ 4-K  , are 
associated  with  the  1st  eigenvalues  At  and  have  nothing  to  do  with  the  2nd  eigenvalues  A2  . From  the  corresponding 
buckled  mode  shapes  ^Yr (£)  shown  in  Figure  8(a)  one  sees  that  the  1st  one  ^Y \(E,)  and  the  3ld  one  ^Y3(^) . take  the 
same  form  as  the  first  and  second  buckled  mode  shapes  of  the  classical  axial-loaded  “1-span”  C-C  beam  with  N = mfV2 , 
cc-Z|  „(£)  and  cc Z2  a ( b ) - respectively,  shown  in  Figure  A. 1(b)  (in  Appendix).  Since  the  “equivalent”  critical  flow 
velocities  corresponding  to  CCZ  lcr{§)  and  ccZlcr{£)  are  given  by  Eq.  (A.10):  ccVlcr  =2  K and  ccV7cr  =4  K , the  1st  and 
3ld  critical  flow  velocities  of  the  pipe  are  given  by  cc^L  cr  ~ 2;T=ccV1*cr  and  cent  er  ~ 4 lt=ccV* cr  as  shown  in  Table  3(a) 
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and  Figure  7(a-l).  From  the  last  table  and  figure,  one  sees  that,  the  pipe  will  be  “buckled”  if  the  flow  velocity  ^'U ] 
(associated  with  the  1st  complex  eigenvalues  A^ ) is  greater  than  the  1st  critical  flow  velocity  (cc^i>cc^iacr)  or  3ld  one 
(cc^i  >cc^wr)»  anc*  wiU  restore  stability  and  coupled  with  the  oscillation  (or  flutter)  associated  with  the  2nd  complex 
eigeanvalyes  A2  if  ccU i is  grater  than  the  2nd  critical  flow  velocity  and  less  than  the  3rd  one  ( cc^i*f>cr<cc^i*<cc^wr  )• 

For  the  classical  axial-loaded  “1-span”  C-C  beam , its  rth  critical  buckled  mode  shape  is  determined  by  Eq.  (A. 11), 
i.e.  ccZr cr(x)  = Ar(l  — cosarx)  . Since  I cos arx  I < 1 .0 , the  magnitude  of  the  last  mode  shape  is  greater  than  or  equal  to 

zero  (i.e.,  ccZrcr(x)>  0)  as  one  may  see  from  Figure  A. 1(b).  This  is  the  reason  why  the  mode  shape  {qCY2  (£) 

corresponding  to  the  2nd  critical  flow  velocity  cc^ibcr  =8.9868  does  not  appear  in  Figure  A.  1(b)  for  the  classical  axial- 
loaded  “1-span”  C-C  beam.  It  is  noted  that  the  Argand  diagrams  shown  in  Figure  7 are  based  on  Eq.  (17),  in  which,  the 
overall  effective  damping  matrix  [c]  is  obtained  from  the  assembly  of  the  element  damping  matrix  [c]f  given  by  Eqs.  (13) 

and  (15a)  with  Coriolis  force  as  the  matrix  coefficients  for  the  case  of  E = Ca=  0.  Thus,  the  mode  shape  ^Y2  (q) 
corresponding  to  the  2nd  critical  flow  velocity  cc^ibcr  =8.9868  is  due  to  the  Coriolis  effect  and  this  effect  should  be  the 
mechanism  for  the  pipe  restoring  stability  and  coupling  with  the  oscillation  (or  flutter)  associated  with  the  2nd  complex 
eigeanvalyes  A2  if  ccuw,cr<ccul<cculc.cr  ■ 

For  the  2-span  C-C  pipe  shown  in  Figure  6(b),  the  Argand  diagrams  associated  with  its  1st  and  2nd  complex 
eigenvalues  are  shown  in  Figure  7(b-l)  and  (b-2),  respectively.  From  Figure  7(b-l)  one  sees  that  the  1st  critical  flow 

velocity  is  given  by  cc^iC  = 8.9869  and  the  corresponding  “buckled”  mode  shape  (denoted  by  the  curve  — • — ) is  shown 
in  Figure  8(b),  which  takes  the  same  form  for  the  2nd  buckled  mode  shape  (cf.  the  curve  — + — ) of  the  1-span  C-C  pipe 
shown  in  Figure  8(a),  i.e.,  (^')~^Y2  {£,) , thus,  ccU*lcr  = 8.9869  ~ cc^ibcr^  as  they  should  be.  Furthermore,  from 

Figure  7(b-2)  one  sees  that  the  2nd  critical  flow  velocity  of  the  2-span  C-C  pipe  is  given  by  $?cU2cr=  An  with 
corresponding  buckled  mode  shape  shown  in  Figure  8(b).  Because  the  last  buckled  mode  shape  (q)  takes  the  same 
form  as  the  2nd  mode  shape  ccZ,(<f)  of  the  classical  axial-loaded  “1-span”  C-C  beam  with  ccV2cr=^ft  as  shown  in 
Figure  A.  1 (b),  thus,  ^U*2  cr  = 4#  =ccV2cr . 

For  the  3-span  C-C  pipe  shown  in  Figure  6(c),  the  Argand  diagrams  associated  with  its  1st  and  2nd  complex 
eigenvalues  are  shown  in  Figs.  7(c-l)  and  (c-2),  respectively,  and  the  coordinates  (AlR  ,AW)  and  (A2R,A2I)  for  some 

relevant  points  on  the  loci  of  the  1st  and  2nd  eigenvalues  are  listed  in  Table  3(b).  From  Table  3(b)  and  Figure  7(c-l)  and 
(c-2)  one  sees  that  the  values  of  Aj ; and  A2 ; decrease  with  the  increase  of  flow  velocity  ^U*r , and  Au  = 0 (but  Z, , = 

105.12)  when  = 1 1.5701 , which  is  the  1st  critical  buckling  flow  velocity  [a  „ associated  with  the  1st  eigenvalues. 
For  the  flow  velocities  greater  than  cc^Lcr  =11-5701,  Au  = 0 and  the  value  of  A, , continuously  decreases  with  the 
increase  of  > an£l  both  Al , and  A 9 , are  equal  to  zero  (Au  = A2I  = 0 ) when  7*  = 15.4087  , which  is  the  1st  critical 
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buckling  flow  velocity  cc^iacr  associated  with  the  2nd  eigenvales.  For  the  flow  velocities  ,^'U]  to  be  greater  than 
ccU*-,acr=  15.4087,  Au  = A2 , =0  until  ^ £7*  =17.3924,  which  is  the  2nd  critical  flow  velocity  associated  with  1st  and  2nd 

complex  eigenvalues,  i.e.,  ^£7 *bcr  = 17.3924=  cc^lbcr  ■ Both  the  F'-mode  and  2nd-mode  vibrations  begin  to  couple  when 
the  flow  velocity  is  greater  than  the  last  critical  value.  After  that,  the  pipe  restores  stability  with  Au  = A2 , = 0.073  and 

AlR  =87.435  = -A1R  , and  Alr  =0  (but  A2I  = 15.968)  when  ® Ur  =18.8496,  which  is  the  3rd  critical  buckling  flow 
velocity  <^£7^ cr  associated  with  the  1 5,1  eigenvalues.  For  the  flow  velocities  to  be  greater  than  cc^lc cr  = 1 8.8496 , Alr  = 0 
and  An , continuously  decreases  with  the  increase  of  c('U‘. , and  both  At , and  An , are  equal  to  zero  (i.e.,  Au  = A2 , = 0 ) 
when  ccU*r  =21.7891,  which  is  the  3rd  critical  buckling  flow  velocity  (x;U2c„  associated  with  the  2nd  eigenvalues.  For 
comparison,  the  buckled  mode  shapes  corresponding  to  the  two  critical  flow  velocities  ( cc^ucr  = H-5701  and 

ccU*bcr  = 17.3924  ) associated  with  the  1st  eigenvales  Al  and  those  ( cc^ia  cr  = 15.4087  and  cc^*ibcr  = 17-3924  ) 
associated  with  the  2nd  eigenvales  A2  are  plotted  in  Figure  8(c).  It  is  seen  that  the  mode  displacements  of  all  the  four 
buckled  mode  shapes  are  in  the  range  from  +1.0  to  -1.0,  i.e.,  — 1.0  < . ( q)  < +1.0 , instead  of  greater  than  or  equal  to  zero 

for  the  buckled  mode  shapes  of  the  classical  axial-loaded  “1-span”  C-C  beam  ( ccZr(x)>  0)  shown  in  Figure  A.l(b).  For 
the  last  reason,  the  foregoing  critical  buckling  flow  velocities  of  the  3-span  C-C  pipe  are  not  equal  to  2m 
(i.e.,  ccU*cr  7=-  2r7T ),  with  r = 1,2,3,... . 

For  the  4-span  C-C  pipe  shown  in  Figure  6(d),  from  the  Argand  diagrams  associated  with  its  1st  and  2nd 
eigenvalues  shown  in  Figure  7(d-l)  and  (d-2),  one  sees  that  the  pipe  will  be  buckled  when  ^£7*  = 14.2995 , which  is  the  1st 

critical  buckling  flow  velocity  eerier  associated  with  the  1st  eigenvalues  At.  For  the  flow  velocities  to  be  greater  than 
cc^lcr  =14.2995,  AiI  is  always  equal  to  zero  (Au  = 0)  and  An , decreases  with  the  increase  of  ^*£7* , and  71, ; =0 
when  ccU*  — 17.9737  , which  is  the  2nd  critical  buckling  flow  velocity  associated  with  the  2nd  eigenvalues,  i.e., 
cc^iacr  =17.9737.  From  the  Argand  diagram  shown  in  Figure  7(d-2),  it  is  seen  that  the  pipe  will  restore  stability  at 

cc^Jibcr = 23.5108  and  be  buckled  again  at  cc^lccr  = 25.1330- Sir  . It  is  different  from  the  first-and-second-mode- 
coupled  vibrations  of  the  3-span  C-C  pipe  (cf.  Figure  7(c-l)  and  (c-2))  that  the  vibration  of  the  4-span  C-C  pipe  is 
associated  with  the  2nd  eigenvalues  only  “after  it  restores  stability”  at  cc^2bcr  ~ 23.5108  . The  1st,  2nd  and  4th  buckled  mode 

shapes  of  the  4-span  C-C  pipe  are  shown  in  Figure  8(d).  It  is  seen  that  the  4th  buckled  mode  shape  of  the  4-span  C-C  pipe 
takes  the  same  form  as  the  4th  one  of  the  classical  axial-loaded  “1-span”  C-C  beam  shown  in  Figure  A.  1(b),  thus,  the 

corresponding  critical  flow  velocity  cc^iccr  is  equal  to  ccV^cr,  i.e.,  cc^lccr  ~^=cc^4cr 
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Figure  6:  The  Four  Uniform  Fluid-Conveying  C-C  Pipes  Subjected  to  Axial  load  N = mfU 2 with  Equal  Span 
Lengths  il=  i 2=  £3=  £ 4 = £0.  The  Other  Legends  are  the  Same  as  Those  of  Ligure  2. 


Table  3:  The  Coordinates  ( Ar  R , Xr  J ) for  Some  Points  on  the  Loci  of  the  1st  and  2nd  Eigen  Values  of  the  C-C  Pipes 

(a)  Lor  1-span  C-C  pipe 


mu" 

CC  W r 

0.05 

... 

6.2831 

6.2832 

mU‘ 

CC  V 1 a,cr 

... 

8.9868 

WU' 

CC  1 b,cr 

8.9869 

9.2 
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... 

12.2 

12.3 

12.5664 

12.5665 

mU * 

CC  ^ 1 c,cr 
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15.0 

4 

A .R 

0 

0 

0 

0.067 
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0 

0 

1.408 

7.422 

0 

0 
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1.556 

K, 

22.372 

0.098 

0 

0 

0 
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10.567 

17.095 

27.173 

19.335 

0.068 

0 

0 

0 

K 

^2  ,R 

0 

0 

0 

0 

0 

0 

0 

0 
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17.646 

17.649 
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^2.1 

61.672 

47.496 

47.495 

26.982 

26.980 

22.682 

17.095 

27.173 

40.149 

43.518 

43.519 

56.355 

(b)  Lor  3-span  C-C  pipe 


i3)u* 

CC  r 

0.05 

11.5700 

11.5701 

mU” 

CC  ^ la.cr 

15.4086 

15.4087 
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17.3924 

# 

17.3925 
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ccK,r 
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0 
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CC  ^ 2c,cr 
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K 
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0 

0 

0 
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0 

0.604 

9.482 
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K 

113.83 
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0 

0 

0 

0 

0 

0 
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0.42 

0 

0 

0 

0 

0 

0 
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0 
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0 

0 

0 
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0 
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K, 

166.22 
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0.545 

0 

0 

0 

0.073 

15.967 

15.968 

2.874 

0 

0 
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Figure  7:  The  Argand  diagrams  for  the  1st  and  2nd  complex  eigenvalues  of  (a)  1-,  (b)  2-,  (c)  3-  and  (d) 
4-  span  fluid-conveying  C-C  pipes  shown  in  Figs.  6(a)-(d)  with  mass  ratio  mf  = 0.5 . In  other 

words,  (a-1)  and  (a-2)  are  those  of  1-span  pipe;  (b-1)  and  (b-2)  are  those  of  2-span  pipe; 

(c-1)  and  (c-2)  are  those  of  3-span  pipe;  and  (d-1)  and  (d-2)  are  those  of  4-span  pipe. 
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(c-1)  (c-2) 


(d-1)  (d-2) 


Figure  7:  (Continued) 
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Istto  3rd  critical  buckled  mode  shapes  of 
tie  C-C  pipe  associated  with  1st  and  2nd 
co  m p le  x e ig  e n val  ue  s 

— 0 — 1st  critical  buckled  mode  shape 

2nd  critical  buckled  mode  shape 

3rd  critical  buckled  mode  shape 

3rd  critical  buckled  mode  shape 
(regaining  stab irty for  3-span  p^e) 


1st.  2nd  and  4th  critical  buckled  mode  shapes 
of4-span  C-C  pipe  corresponding  to 
1st  and  2nd  complex  eigen  values 

— 9 — 1st  critical  buckled  mode  shape 

— ^ — 2nd  critical  buckled  mode  shape 

— 0 — 4th  critical  buckled  mode  shape 


Figure  8 The  1st  to  4th  “buckled”  Mode  Shapes  of  (a)  1-,  (b)  2-,  (c)  3-  and  (d)  4-  span 

Fluid-Conveying  C-C  pipes  Shown  in  Figure  6(a)-(d)  with  Mass  Ratio  m*f  = 0.5 . 


3.2.3  With  C-P  end  supports  (orBCs) 

The  configurations  of  the  1-  to  4-  span  fluid-conveying  C-P  pipes  with  equal  span  lengths  l o = 6m  are  shown  in 
Figure  9 and  the  Argand  diagrams  associated  with  the  1st  and  2nd  complex  eigenvalues  of  the  g-span  pipes  are  shown  in 
Figure  1 0(^-1)  and  (q- 2),  respectively,  with  span  number  q = a = 1 , q = b = 2 , q = c = 3 and  q = d = 4 . The  coordinates 
( Ar  R , Ar , ) for  some  relevant  points  on  the  loci  of  the  rth  eigenvalues  (with  r = 1,2 ) for  the  1 - and  2-  span  C-P  pipes  are 
shown  in  Tables  4(a)  and  (b),  respectively.  It  is  similar  to  the  1-span  P-P  pipe  (Figure  4(a))  and  the  1-span  C-C  pipe  (Figure 
7(a))  that  all  the  lowest  three  critical  flow  velocities  of  the  C-P  pipe,  cr^hc,  =4.4935  , a>U]b  cr  = 7.7252  and 
cpU iccr  =10.9042,  are  associated  with  the  1st  eigenvalues  and  have  nothing  to  do  with  the  2nd  ones  as  one  may  see  from 
Table  4(a)  and  Figure  10(a-l)  and  (a-2).  The  pipe  will  be  buckled  if  the  flow  velocity  '('pU * is  greater  than  the  1st  critical 
one  ( cpU Ul  cr)  or  3ld  critical  one  ( cpUlccr  X and  will  regain  stability  if  the  value  of  ppU].  is  greater  than  the  2nd  critical  flow 
velocity  (a>U*lhcr)-  In  addition,  the  vibration  of  the  pipe  is  first-and-second-mode-coupled  for  the  flow  velocity  * to 
be  in  the  range  of  2nd  to  3rd  critical  ones,  i.e.,  cpUlbcr  ^ cp^i  — cp^iccr-  Since  the  corresponding  three  buckled  mode 
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shapes  shown  in  Figure  11(a)  take  the  same  forms  as  the  lowest  three  ones  of  the  classical  axial-loaded  “1-span”  C-P  beam 
shown  in  Figure  A.  1(c),  one  has:  ^.U*Ucr  = 4.4935 =CPV*cr , cpUlb,cr  = 7 -7 252= cpKcr  > cpUu,cr  = 10.9042=CPV3*(T  ■ 

For  the  2-span  C-P  pipe,  the  coordinates  ( Ar  R , Ar , ) for  some  relevant  points  on  the  locus  of  the  rth  complex 

eigenvalue  Ar  show  in  Table  4(b)  and  the  Argand  diagrams  shown  in  Figure  10(b-l)  and  (b-2)  reveal  that  the  pipe  will  be 

buckled  in  1st  mode  if  the  flow  velocity  ^t/*  is  greater  than  the  1st  critical  one  cPU‘la a =7.1498  associated  with  the  1st 

eigenvalue  At  and  in  2nd  mode  if  the  value  of  is  greater  than  the  2nd  critical  one  (cPU2acr  =10.8852  associated  with 

the  2nd  eigenvalue  A2 . The  last  two  buckled  modes  will  restore  stability  when  cpU*.  >11.7554,  thus,  11.7554  is  the  3rd 

critical  flow  velocity  associated  with  both  the  1st  and  2nd  eigenvalues,  i.e.,  cv* ^ u>  cr = cr  ^ 2b  cr  =11-7554.  For  the  flow 

velocities  to  be  greater  than  the  last  critical  values,  the  1st  mode  vibration  and  the  2nd  one  will  be  coupled.  Since  the  1st  and 
2nd  buckled  mode  shapes  shown  in  Figure  11(b)  look  like  the  2nd  and  3ld  buckled  mode  shapes  of  the  classical  axial-loaded 
“1-span”  C-P  beam  shown  in  Figure  A.  1(c)  to  some  degree,  respectively,  the  corresponding  critical  flow  velocities 

a’^uicr  =7.1498  and  a>U*2aa  = 10.8852  are  somewhat  close  to  CIV2cr  - 7.7253  and  CPV2cr  = 10.9041,  respectively. 

For  the  3-span  C-P  pipe,  from  its  Argand  diagrams  shown  in  Figure  10(c-l)  and  (c-2)  one  sees  that  the  pipe  will 
be  buckled  in  1st  mode  if  the  flow  velocity  ^t/*  is  greater  than  the  1st  critical  one  cpU\cr  =10.0379  associated  with  the 

1st  eigenvalue  A , and  in  2nd  mode  if  the  value  of  cpU*  is  greater  than  the  2nd  critical  one  cPU2cr  = 13.4803  associated  with 
the  2nd  eigenvalue  A 0 . From  the  lowest  two  buckled  mode  shapes  shown  in  Figure  11(c)  one  sees  that  the  1st  one  look  like 
the  3rd  buckled  mode  shapes  for  the  classical  axial-loaded  “1-span”  C-P  beam  shown  in  Figure  A. 1(c)  to  some  degree,  thus, 
the  corresponding  critical  flow  velocity  cp^icr  = 10.0379  is  somewhat  close  to  CPV2cr  = 10.9041. 

For  the  4-span  C-P  pipe  with  its  Argand  diagrams  shown  in  Figure  10(d-l)  and  (d-2),  it  is  seen  that  the  lst-mode 
buckling  occurs  when  the  flow  velocity  is  greater  than  the  1st  critical  one  cpU*lcr  =13.0366  associated  with  the  1st 

eigenvalue  A{  and  the  2nd-mode  buckling  occurs  when  the  value  of  ^PU*r  is  greater  than  the  2nd  critical  one 
cpU*2cr  =16.0408  associated  with  the  2nd  eigenvalue  A2 . Since  the  corresponding  lowest  two  buckled  mode  shapes  shown 
in  Figure  11(d)  are  different  from  anyone  of  the  lowest  four  buckled  mode  shapes  for  the  classical  axial-loaded  “1-span” 
C-P  beam  shown  in  Figure  A.  1(c),  the  foregoing  two  critical  flow  velocities,  cpU\cr  =13.0366  and  cPU2 cr  =16.0408,  are 

not  close  to  any  of  the  equivalent  ones,  CPVrcr  ( r = 1 — 4),  given  by  Eqs.  (A.  16)  and  (A.  19).  Flowever,  since  the  1st  critical 
buckled  mode  shape  of  the  4-span  C-P  pipe  shown  in  Figure  11(d),  ^Tlcr,  looks  like  the  4th  buckled  mode  shape  of  the 
classical  axial-loaded  “1-span”  P-P  beam  shown  in  Figure  A.  1(a),  PPZ4cr,  its  critical  flow  velocity  cfoicr  =13.0366  is 

slightly  greater  than  PPV4cr  = \7t  = 12.56  , because  the  stiffness  of  the  C-P  pipe  is  slightly  greater  than  that  of  the 
corresponding  P-P  pipe. 
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Figure  9:  The  Four  Uniform  Fluid-Conveying  C-P  Pipes  Subjected  to  Axial  Load  N = mfU2  with  Equal 
Span  Lengths  £l=£2=£3=£A=£0.  The  Other  Legends  are  the  same  as  those  of  Ligure  2 


Table  4:  The  Coordinates  (ArR,  Xr , ) for  Some  Relevant  Points  on  the  Loci  of  the  1st  and  2nd  Complex 
Eigen  Values  of  the  Lluid-Conveying  C-P  Pipes:  (a)  1-Span;  (b)  2-Span 
(a)  1-Span  C-P  Pipe 


mu' 

CP^  r 

0.05 

4.4934 

4.4935 

mU’ 

Cl’u  ill  .IT 

7.7252 

"V 

CP  ^ lb.i 

7.7253 

7.9 

8.0 

10.7 

10.8 

10.9041 

10.9042 

CP  ^ Ic.cr 

12.0 

A 

A.X 

0 

0 

0 

0.102 

0.159 

0 

-7.598 

9.655 

15.887 

12.637 

0 

0.241 

1.802 

A, 

15.417 

0.029 

0 

0 

0 

0.197 

11.170 

11.756 

15.162 

14.966 

0.397 

0 

0 

0 

A 

^ l.R 

0 

0 

0 

0 

0 

0 

0 

7.598 

-9.655 

-15.887 

-12.637 

0 

0 

27.367 

A.l 

49.964 

41.655 

41.655 

14.597 

14.593 

11.170 

11.756 

15.162 

14.966 

23.543 

23.574 

36.829 

(b)  2-Span  C-P  Pipe 


mu * 

CP^  r 

0.05 

... 

7.1497 

7.1498 

'AK.cr 

7.3 

... 

10.8851 

10.8852 

CpKa.cr 

... 

11.7554 

# 

11.7555 

11.9 

12.0 

A 

A ,R 

0 

0 

0 

0.291 

2.023 

4.354 

4.354 

4.108 

36.311 

-37.530 

-38.299 

A. , 

46.055 

0.135 

0 

0 

0 

0 

0 

0 

0.123 

4.712 

6.031 

A 

^2  ,R 

0 

0 

0 

0 

0 

0 

0 

-0.307 

-4.105 

-36.311 

37.530 

38.299 

A, 

79.684 

58.255 

58.254 

57.233 

0.255 

0 

0 

0.125 

4.712 

6.031 
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Figure  10:  The  Argand  Diagrams  for  the  1st  and  2nd  Complex  Eigen  Values  of  (a)  1-,  (b)  2-,  (c)  3-  and  (d)  4-  Span 
Fluid-Conveying  C-P  Pipes  Shown  in  Figs.9(a)-(d)  with  Mass  Ratio  nif  = 0.5  . In  Other  Words,  (a-1) 

and  (a-2)  are  Those  of  1-span  pipe;  (b-1)  and  (b-2)  are  those  of  2-span  pipe;  (c-1)  and  (c-2)  are  those  of 
3-Span  Pipe;  and  (d-1)  and  (d-2)  are  those  of  4-Span  Pipe 
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Figure  10:  (Continued) 
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Figure  11:  The  1st  to  3rd  Buckled  Mode  Shapes  of  (a)  1-,  (b)  2-,  (c)  3-  and  (d)  4-  Span 

Fluid-Conveying  C-P  Pipes  Shown  in  Figures.  9(a)-(d)  with  Mass  Ratio  mf  = 0.5 

3.3  COMPARISON  BETWEEN  1st  CRITICAL  FLOW  VELOCITIES  FOR  P-P,  C-C  AND  C-P  PIPES 

For  comparisons,  the  1 M critical  flow  velocities  for  the  1-,  2-,  3-  and  4-  span  fluid-conveying  pipes  with  P-P,  C-C 
and  C-P  end  conditions  are  listed  in  Table  5.  In  which,  the  1st  critical  flow  velocities  shown  in  the  3rd  row  for  the  P-P  pipe, 

W-Ulcr  (with  q = 1,2, 3, 4),  are  obtained  from  the  Argand  diagrams  shown  in  Figure  4(a-l),  (b-1),  (c-1)  and  fd-1), 
respectively;  those  shown  in  the  4th  row  for  the  C-C  pipe,  eerier  (with  q = 1, 2,3,4),  are  obtained  from  Figure  7(a-l),  (b-1), 

(c-1)  and  fd-1),  respectively;  and  those  shown  in  the  5th  row  for  the  C-P  pipe,  cpUlcr  (with  q = 1,2, 3, 4),  are  obtained  from 
Figure  10(a-l),  (b-1),  (c-1)  and  (d-1),  respectively.  From  Table  5 one  sees  that: 

• For  a ‘'(/-span”  P-P  pipe  with  equal  span  lengths,  its  1st  critical  flow  velocity  of  buckling  is  given  by  ppU^cr  = q7t , 
which  is  equal  to  the  rth  equivalent  critical  flow  velocity  obtained  from  the  theory  of  static  stability 
(for  the  classical  axial-loaded  “1-span”  P-P  beam),  PPV*cr  = rn  , i.e.,  ppU^ cr= PPV* cr  = m (with  q = r).  Since  the 
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expression  PPVrcr  = rn  obtained  from  the  theory  of  static  stability  is  independent  on  the  mass  ratio  mf,  so  is 

the  expression  ■ „pUlcr=PPVrcr  = t'7l  (with  q = r).  The  last  relationship  is  not  available  for  the  “(/-span”  (q>  1) 
pipe  with  the  other  end  conditions  (such  as  C-C  or  C-P). 

• For  an  “1-span”  pipe  with  S end  conditions  (with  .V  = P-P,  C-C  or  C-P),  its  1st  critical  flow  velocity  obtained  from 
the  theory  of  dynamic  stability  presented  in  this  paper,  gUlcr , is  equal  to  the  1st  equivalent  one  obtained  from  the 

theory  of  static  stability  shown  in  the  Appendix  (of  this  paper),  sVlcr , i.e.,  sUicr=sVlcr . 

• Among  the  1st  critical  flow  velocities  of  the  “1-span”  P-P,  C-C  and  C-P  pipes,  the  order  of  the  magnitudes  is 
i'pU \ cr<cpU \ ,-r<ccU u.r , this  is  reasonable  because  the  stiffness  of  the  P-P  pipe  is  smallest  and  that  of  the  C-C 

pipe  is  largest.  The  last  trend  is  also  true  for  the  “g-span”  pipes,  i.e.  // U *cr < // U *cr < //? U *c:r  with  q > 2 , and  the 
divergences  between  them  decrease  with  the  increase  of  span  number  q. 


Table  5:  The  “first”  Critical  Dimensionless  flow  Velocities  sUlcr  of  the  “(/-span”  Fluid-Conveying 

Pipe  with  S BCs  and  the  Lowest  “four”  Equivalent  Critical  ones  sV*cr  ( /-  = 1 ~ 4 ) of  the 
Classical  Axial-Loaded  “1-span”  beam  with  S BCs 


BCs  for  r/-Span  Pipe 
S 

“First  “ Critical  Dimensionless  Flow  Velocities  of  “(/-Span”  Pipe,  ''(T/I  r 

{VKr,Cq  = \) 

™Ulr(q  = 2) 

°1Km= 3) 

= 4) 

P-P 

71 

2 n 

3 n 

An 

C-C 

In 

8.9869 

11.5701 

14.2995 

C-P 

4.4935 

7.1498 

10.0379 

13.0366 

BCs  for  axial-loaded 
1-span  ( q = 1 ) beam 

S 

Lowest  “four”  equivalent  critical  dimensionless  flow  velocities  of  “1-span”  beam , sV*cr 

sVlUbr  = 1) 

sKcr(r  = 2') 

sVir(r  = 3) 

sV;,cr(r  = A) 

P-P 

7T  (=3.1416) 

2 ?r  (=6.2832) 

3 n (=9.4248) 

4 ^(=12.5664) 

C-C 

2 ;r  (=6.2832) 

4 ;r  (=12.5664) 

6 ;r  (=18.8496) 

8 7T  (=25. 1 328) 

C-P 

4.4935 

7.7253 

10.9041 

14.0662 

a q = total  number  of  spans  for  each  pipe 


b r = mode  number  for  a specified  vibration  of  beam 

3.4  FREE  VIBRATION  ANALYSES  OF  A 4-SPAN  PIPE  WITH  VARIOUS  END  CONDITIONS 

In  general,  the  1st  critical  flow  velocities  obtained  in  the  last  subsection  are  much  higher  than  the  realistic  ones, 
thus,  this  subsection  illustrates  a practical  example  to  show  the  influence  of  fluid  pressure  intensity  p f , flow  velocity  U 

and  Coriolis  force  on  the  lowest  four  natural  frequencies  0Jr  of  a 4-span  fluid-conveying  pipe  with  P-P,  C-C  and  C-P  end 
conditions  as  shown  in  Figure  2(d),  6(d)  and  9(d),  respectively.  The  dimensions  and  physical  properties  of  the  pipes  are  the 
same  as  those  mentioned  at  the  beginning  of  this  section:  outer  diameter  da  = 0.254  m,  inner  diameter  dt  = 0.2286  m, 

Young’s  modulus  £ = 2.068xl0n  N/m2  , mass  density  of  pipe  pp  =7850  kg/m3  , mass  density  of  water  (outside  the 
pipe)  pw  =1000  kg/m3  and  total  length  of  the  pipe  L = 6x4  = 24  m.  Besides,  we  assume  that  the  mass  density  of  fluid 
(inside  the  pipe)  is  pf  =1000  kg/m3  and  the  externally  applied  axial-load  is  Pp  =1.0xl06N.  Three  cases  are  studied: 
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Case  1 is  with  p f = U =0  (no  fluid  pressure  and  no  Coriolis  force);  Case  2 is  with  p , = 1.0x10s  N/m2  and  U = 20  m/s 
(with  Coriolis  force  neglected);  Case  3 is  the  same  as  Case  2 except  that  the  Coriolis  force  is  considered.  The  results  are 
listed  in  Table  6 with  total  elements  for  each  pipe  to  be  ne  =50x4  = 200  . Since  each  node  has  two  DOFs,  the  total 
effective  DOFs  for  the  4-span  P-P,  C-C  and  C-P  pipes  are  397,  395  and  396,  respectively.  Among  the  three  pipes,  the  total 
effective  DOFs  of  the  P-P  pipe  is  largest  so  that  its  stiffness  is  smallest,  on  the  contrary,  the  total  effective  DOFs  of  the  C-C 
pipe  is  smallest  so  that  its  stiffness  is  largest.  This  is  the  reason  why  the  order  for  the  lowest  four  natural  frequencies  of  the 
three  pipes  is:  %!(Or  < // CO r < a'0Jr  (with  r = 1-4)  as  one  may  see  from  Table  6 for  either  Case  1,  Case  2 or  Case  3.  From 

Table  6 one  also  sees  that  either  the  “effective  axial  load”  due  to  fluid  pressure  pfAf  as  well  as  centrifugal  force  mfU 2 or 

the  “Coriolis  force’’  can  reduce  the  lowest  four  natural  frequencies  of  the  fluid-conveying  pipe,  however,  the  effect  of  the 
“Coriolis  force”  is  much  smaller  than  that  of  the  “effective  axial  load”. 

Table  6:  Influence  of  Fluid  Pressure  Intensity  pf  , Flow  Velocity  V and  Coriolis  Force  on  the  Lowest  Four 
Natural  Frequencies  cor  of  a 4-Span  Fluid-Conveying  Pipe  with  P-P,  C-C  and  C-P  End  Conditions 
Shown  in  Figure  2(d),  6(d)  and  9(d).  ( pf  =1000  kg/m3  and  Pp  =1.0xl06N) 


Case 

Pressure 

Intensity 

Pf  ( N/m2  ) 

Flow 
Velocity 
U (m/s) 

Is  Coriolis 
force 

considered? 

End 

conditions 

S 

Natural  Frequencies  0)r  (Rad/Sec) 

co, 

co2 

co 

CO 4 

1 

(Undamped) 

0 

0 

No 

P-P 

69.9263 

84.6625 

118.3916 

156.5799 

C-C 

84.6625 

118.3916 

156.5799 

177.4897 

C-P 

73.8766 

100.1403 

137.8284 

171.4725 

2 

(Undamped) 

1.0x10s 

20 

No 

P-P 

69.6854 

84.4540 

118.2250 

156.4430 

C-C 

84.4540 

118.2250 

156.4430 

111  3112 

C-P 

73.6458 

99.9546 

137.6775 

171.3480 

3 

(Damped) 

1.0x10s 

20 

Yes 

P-P 

69.6782 

84.4457 

118.2152 

156.4319 

C-C 

84.4457 

118.2152 

156.4318 

177.3588 

C-P 

73.6383 

99.9455 

137.6670 

171.3361 

Figure  12(a)-(c)  show  the  lowest  four  natural  mode  shapes  of  the  4-span  pipe,  <4/yr(x)  ( r = 1 — 4).  In  which. 
Figure  12(a)  is  for  the  P-P  pipe,  (ppYr(x)  • Figure  12(b)  is  for  the  C-C  pipe,  //T,. (x) , and  Figure  12(c)  is  for  the  C-P  pipe, 

cpK(x)  ■ It  is  noted  that  the  solid  curves  ( ) with  symbols,  O,  X,  A and  □,  denote  the  1st,  2nd’  3ldand  4th  un-damped 

mode  shapes  (for  Case  1),  respectively,  while  the  dashed  curves  ( ) with  symbols,  •,  +,  ▲ and  ■,  denote  the 

corresponding  damped  ones  (for  Case  3),  respectively.  Since  the  un-damped  natural  frequencies  are  very  close  to  the 
corresponding  damped  ones  as  shown  in  Table  6,  so  are  the  associated  natural  mode  shapes  as  one  may  see  from  Figure  12. 
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(a) 


The  lowest  4 natural  mode  shapes  of  the  4-span  pipe 
— © — 1st  mode  shape  (un-damped:  Case  1)  — • — - 1st  mode  shape  (damped:  Case  3) 

— X — 2nd  mode  shape  (un-damped:  Case  1)  — - 2nd  mode  shape  (damped:  Case  3) 

— — 3rd  mode  shape  (un-damped:  Case  1)  — — 3rd  mode  shape  (damped:  Case  3) 

— O — 4th  mode  shape  (un-damped:  Case  1)  — - 4th  mode  shape  (damped:  Case  3) 


10  11  12  13  14  15  16  17  18  19  20  21  22  23  24 
Axial  coordinates,  x(m) 


(b) 


Axial  coordinates,  x(m) 


Figure  12:  The  Lowest  Four  Natural  Mode  Shapes  of  a 4-Span  Pipe  with  Three  End  Conditions:  (a)  P-P, 

(b)  C-C  and  (c)  C-P.  Where  The  solid  Curves,  — O — , — X — , — A — and  — □ — , denote  the 
1st,  2nd'  3rd  and  4th  un-Damped  Mode  Shapes  (Case  1),  while  the  Dashed  Curves,  — • — , — +— , 

— A—  and  — ■ — , Denote  the  Corresponding  Damped  Ones  (Case  3). 

4.  CONCLUSIONS 

Based  on  the  foregoing  numerical  results  and  discussions,  we  draw  the  following  conclusions: 

• For  an  “1-span”  fluid-conveying  pipe  with  S end  conditions  (with  S = P-P,  C-C  or  C-P),  its  “1st”  critical  flow 
velocity  obtained  from  the  theory  of  dynamic  stability  presented  in  this  paper,  l !?f/l*(r , is  equal  to  the  “1st” 
equivalent  one  obtained  from  the  theory  of  static  stability  (for  the  classical  axial-loaded  “1-span”  beam ) shown  in 
the  Appendix,  sVlcr,  i.e.,  ^sUlcr=sVlcr . This  is  a reasonable  result,  because,  according  to  the  theory  of  either 
dynamic  or  static  stability,  the  first  critical  flow  velocity  is  obtained  from  the  condition  that  the  first  natural 
frequency  of  the  vibrating  system  approaches  zero  (i.e.  ty,  = 0 ).  This  reasonable  result  may  also  be  one  of  the 
evidences  that  the  theories  presented  in  the  text  and  the  appendix  of  this  paper  should  be  reliable. 
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• Since  the  “1st”  critical  buckled  mode  shapes  of  the  one-,  two-,  three-  and  four-  span  fluid -conveying  “P-P  pipes ” 
with  identical  span  lengths,  {PPYlcr (£)  (q  = 1-4  ),  take  the  same  forms  as  the  1st,  2nd,  3rd  and  4th  buckled  mode 

shapes  of  the  classical  axial-loaded  “1-span”  P-P  beam,  PPZr  cr(^)  ( r = 1 — 4),  respectively,  and  the  equivalent 

“rth”  dimensionless  critical  flow  velocity  of  buckling  for  the  latter  beam  is  given  by  PPVr  cr  = rJZ  , the 
corresponding  “1st”  dimensionless  critical  flow  velocity  of  buckling  for  the  “(/-span”  P-P  pipe  (with  identical  span 
lengths)  is  given  by  PPU  lcr=PPV*cr  = q7t  (with  q = r).  The  last  relationship  (PPU*Ur  = qn  is  not  available  for  the 
“multi- span”  C-C  or  C-P  pipe. 

• For  the  “1-span”  fluid-conveying  pipe  with  P-P,  C-C  and  C-P  BCs,  the  stiffness  of  the  P-P  pipe  is  smallest  and  that 
of  the  C-C  pipe  is  largest,  so  that  the  magnitude  order  for  the  associated  “1st”  critical  flow  velocities  is: 

ppU\cr<cpU\cr<ccU\cr  • This  conclusion  is  also  true  for  the  “g-span”  fluid-conveying  pipe,  i.e. 

/>/*  U i c,  <(cr  U ur<cx  U icr  with  q>  2,  and  the  divergences  between  them  decrease  with  the  increase  of  the  total 
span  number  q. 

• For  the  “1-span”  P-P,  C-C  or  C-P  fluid-conveying  pipe,  all  three  critical  dimensionless  flow  velocities 

( jt/ lacr,  sUiba-  and  s^iccr)  are  associated  with  the  “1st”  complex  eigenvalues  (A:).  It  is  found  that  the  pipe 
will  be  buckled  (and  lose  stiffness)  if  the  flow  velocity  aRUl  is  greater  than  the  “1st”  or  “3rd  critical  one 
(i.e.,  aX>asU*lacr  or  t/*CCI.)  and  will  restore  stability  (or  regain  stiffness)  if  the  flow  velocity  is  greater 

than  the  2nd  critical  one  (i.e.,  sU1>gUlbcr).  Since  the  pipe  regains  its  stiffness  when  sU1>gUlbcr  and  will  be 

buckled  again  when  gU1  'Ulccr , this  may  be  the  reason  why  the  imaginary  part  of  the  “1st”  complex 
eigenvalue  ( Al , ) is  equal  to  that  of  the  “2nd”  one  (A2I),  i.e.,  Au  = A2 , , and  the  real  part  of  the  “1st”  complex 
eigenvalue  (AlR)  is  “opposite”  to  that  of  the  “2nd"  one  ( A1 R ),  i.e.,  AlR  = —A2  R , if  the  flow  velocity  is 

greater  than  the  2nd  and  less  than  the  3rd  critical  ones  by  about  0.2,  i.e.,  C^^ihcr  +0.2)  < ^ U <({isUlccr  —0.2)  . 
In  other  words,  for  the  last  case,  the  oscillations  of  the  pipe  are  first-and-second-mode-coupled. 

• Except  for  the  3-span  C-C  pipe,  for  a “(/-span”  P-P,  C-C  or  C-P  pipe  with  q > 2 , the  1st  critical  flow  velocity 
((?tO  is  only  associated  with  the  1st  complex  eigenvalue  (At)  and  the  2nd  critical  flow  velocity  ( sU2cr ) is 
only  associated  with  the  2nd  complex  eigenvalue  ( A2 ),  so  that  the  corresponding  critical  buckled  mode  shapes, 
(9s^i,cr(^)  and  (9s*vr(£> . are  not  coupled. 

• The  main  difference  between  the  critical  flow  velocity  leading  to  buckling  of  the  fluid-conveying  pipe  obtained 
from  the  theory  of  dynamic  stability  presented  in  this  paper,  gUrcr,  and  the  “equivalent”  one  obtained  from  the 

theory  of  static  stability  (for  the  classical  axial-loaded  “1-span”  beam ) shown  in  the  Appendix,  sV*cr,  is  the 
Coriolis  force  arising  from  the  flowing  fluid,  thus,  the  buckled  mode  shape  associated  with  the  “1st”  complex 
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eigenvalue  appearing  in  the  fluid-conveying  pipe  and  missing  in  the  classical  axial-loaded  “1-span”  beam  is  due 
to  the  Coriolis  effect,  and  so  is  the  first-and-second-mode -coupled  vibrations  of  the  fluid-conveying  pipe.  In  other 
words,  the  information  about  restoring  stability  of  the  pipe  from  its  buckling  state  and  coupling  oscillations  of  the 
1 sl  and  2nd  modes  can  only  be  obtained  from  the  theory  of  dynamic  stability. 

• For  the  3-span  C-C  pipe,  although  the  pipe  restores  its  stiffness  at  the  second  critical  flow  velocity  ( cc^u,cr ) 
associated  with  the  “1st”  complex  eigenvalue  ( At ) and  the  latter  is  equal  to  the  second  critical  flow  velocity 
( cc^it  cr ) associated  with  the  “2nd”  complex  eigenvalue  ( A2 ),  i.e.  cc^u>  cr  = cc  1 c,-  = 17.3924 , it  is  found  that, 
in  the  last  condition,  the  corresponding  buckled  mode  shape  associated  with  Ax  is  different  from  that  associated 
with  A 2 as  one  may  see  from  Figure  8(c).  This  result  is  reasonable,  because  the  (buckled)  mode  shapes  are 
dependent  on  the  corresponding  eigenvalues  ( Ax  and  A,, ) and  Al  = Al  R + iAx , in  not  equal  to  A2  = A2  R + iA2 , (i.e., 
Al  ± A2 ) as  one  may  see  from  column  11  of  Table  3(b). 

• Besides  the  stability  analysis,  the  presented  FEM  is  also  available  for  studying  the  influence  of  certain  parameters 
on  the  free  vibration  characteristics  of  a multi-span  fluid-conveying  pipe  subjected  to  arbitrary  externally  applied 
axial-load  Pp  , internal  fluid  pressure  intensity  p f and  fluid  velocity  U with  various  end  supports  (or  BCs).  For 

the  4-span  fluid-conveying  pipe  with  P-P,  C-C  and  C-P  BCs  studied,  it  has  been  found  that  either  the  “effective 
axial  load”  due  to  fluid  pressure  as  well  as  centrifugal  force  or  the  “Coriolis  force”  can  reduce  the  lowest  four 
natural  frequencies,  however,  the  effect  of  the  “Coriolis  force”  is  much  smaller  than  that  of  the  “effective  axial 
load”. 

APPENDIX 

Equivalent  Critical  Flow  Velocities  and  Buckled  Mode  Shapes  for  a Single-Span  Fluid-Conveying  Pipe  Obtained 
from  Theory  of  Static  Stability 

A.l  Axial-Foaded  1-Span  P-P  Beam 

According  to  the  theory  of  static  stability  [30],  the  rth  critical  buckling  load  of  an  axial-loaded  1-span  P-P  beam 
(in  air)  is  given  by 

rPNr,cr=EI(m/L)2  (A.l) 

For  the  case  with  P = p fAf  = 0 , from  Eq.  (16c)  one  obtains  the  equivalent  axial  load 

N = mfV2  (A. 2) 

Where  the  flow  velocity  U is  replaced  by  the  equivalent  one  V , for  convenience. 

From  Eqs.  (A.l)  and  (A. 2),  one  obtains  the  rth  equivalent  critical  flow  velocity  of  buckling  for  a single-span  P-P 
pipe  to  be 

PPVr„=(rx/L)jEI/m f (A.3) 

According  to  Eq.  (23b),  the  corresponding  rth  dimensionless  critical  flow  velocity  is  given  by 
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eeVlr  = Vr,crLjmf/(EI ) = rK  (A.4) 

Furthermore,  from  Ref.  [30],  one  obtains  the  following  rth  critical  buckled  mode  shape 

PPZrCr(x)  = Arsinarx  (A. 5) 

Since 

ar=JppNr,cr/(EI)=Y7r/L  (A.6) 

The  substitution  of  Eq.  (A.6)  into  Eq.  (A.5)  produces 

pPZrCr(x)  = Arsin(rm/L)  (A.7) 


For  Ar  = 1.0  , one  has  lowest  four  critical  buckled  “unit-amplitude”  mode  shapes  as  shown  in  Figure  A.  1(a). 


(a) 


(b) 


• 

1st  mode  with  PPV*cr  = it 

— 1 — 

2nd  mode  with  PPV2*cr  = 2ic 

▲ 

3rd  mode  with  pPV3*cr  = 3n 

-H-- 

4th  mode  with  pPV4*,,r  = 4ic 

0.3  0.4  0.5  0.6  0.7 

Axial  coordinates,  x/L 


1st  mode  with  cc  ,cr  - 2jc 


4th  mode  with  CCV4  cr  = 8n 


1st  mode  with  cpV,  cr  = 4.4935 
2nd  mode  with  CPV2)cr  = 7.7253 
3rd  mode  with  CPV3  cr  = 1 0.9041 
4th  mode  with  cpV^.  = 14.0662 


0.4  0.5  0.6 

Axial  coordinates,  x/L 

Figure  A.l:  The  Lowest  4 Critical  Buckled  “Unit-Amplitude”  Mode  Shapes  of  a Uniform  Axial-Loaded  1-Span 
(a)  P-P,  (b)  C-C  and  (c)  C-P  Beams  with  the  Curves  Attached  by  SYMBOLS,  •,  +,  ▲ and  ■, 
Denoting  the  1st,  2nd,  3rd  and  4th  Ones,  Respectively 

A.2  Axial-Loaded  1-Span  C-C  Beam 

The  rth  critical  buckling  load  of  an  axial-loaded  1-span  C-C  beam  is  given  by  [30] 
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ccNr,cr=4EI(rx/L)2  (A.8) 

For  the  case  withP^  = pfAf  =0  , from  Eq.  (16c)  one  obtains (V  = mfV 2 , thus,  the  equivalent  rth  critical  flow 


velocity  of  buckling  is  given  by 

CcVr,cr=2(m/L)jEl/mf  (A.9) 

Thus,  from  Eq.  (23b),  the  corresponding  rth  dimensionless  critical  flow  velocity  is 

ccV:,cr=ccVr,crL^mf/(EI ) = 2r;r  (A.  10) 

Furthermore,  from  [30],  one  obtains  the  rth  critical  buckled  mode  shape  to  be 

ccZrcr(x)  = Ar  (1-  cos  ar  x)  (A.  1 1) 

Since 

ar  = yj  ppNrcr/(EI)  = IrKjh  (A.  12) 

The  substitution  of  Eq.  (A.  12)  into  Eq.  (A.  11)  leads  to 

ccZrcr(x)  = Ar[l-cos(2m/L)]  (A.  13) 


For  Ar  = 0.5  , one  has  lowest  four  critical  buckled  “unit-amplitude”  mode  shapes  as  shown  in  Figure  A.  1(b). 

A.3  Axial-loaded  1-span  C-P  beam 

The  rth  critical  buckling  load  of  an  axial-loaded  1-span  C-P  beam  is  given  by  [33] 

CPNr,cr=£;EI  (A.  14) 

For  the  case  with  Pp  = pfAf  =0,  from  Eq.  (16c)  one  obtains  N = mfV2,  thus,  the  equivalent  rth  critical  flow 
velocity  of  buckling  is  given  by 

crK,cr=C-  4EIlmf  (A.  15) 

Now,  from  the  above  equation  and  Eq.  (23b),  one  obtains  the  corresponding  rth  dimensionless  critical  flow 

velocity 

cpKr=cpVr^mf /(El)  = £L  (A.  16) 

Furthermore,  from  Ref.  [33],  one  obtains  the  following  rth  critical  buckled  mode  shape 

CpZ  rcrix)  = Ar[S\n£  rL%  - £ rLcose  rU;  + £ rL(\  - £)]  (A.  17) 

Where  £ = x/L  . 

Since  Figure  2.17  of  Ref.  [33]  reveals  that 
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SrL  = 4.4935, 7.7253,  10.9041,  14.0662  (for  r = 1-4)  (A.  18) 

from  Eqs.  (A. 17)  and  (A. 18)  one  obtains  the  lowest  four  unit-amplitude  critical  buckled  mode  shapes  of  the  C-P 
beam  as  shown  in  Figure  A.  1(c). 
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